Scaling Function, Universality and Analytical Solutions of Generalized One-Species 

Population Dynamics Models 
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We consider several one-species population dynamics model with finite and infinite carrying ca- 
pacity, time dependent growth and effort rates and solve them analytically. We show that defining 
suitable scaling functions for a given time, one is able to demonstrate that their ratio with respect 
to its initial value is universal. This ratio is independent from the initial condition and from the 
model parameters. Although the effort rate does not break the model universality it produces a 
transition between the species extinction and survival. A general formula is furnished to obtain the 
scaling functions. 



I. INTRODUCTION 

The population growth problem is one of the most im- 
portant of the scientific knowledge. It is crucial not only 
to describe ecological systems like bacterial or virus pop- 
ulations, but also to understand economic behaviors and 
human population growth 

The simplest way to deal with population growth is to 
consider that their individuals do not interact with exter- 
nal ones. This is represented by the so-called one-species 
population dynamics models. These models quantify 
the size (number of individuals) N(t) > of a given 
population, at a certain time t, given its: initial size 
No = N(0) > 0; growth rate k > and the environmen- 
tal carrying capacity K = JV(oo) > 0. The environmen- 
tal carrying capacity takes into account all possible in- 
teractions among individuals, species and resources into 
a single parameter. If one assumes that the population 
lives in an environment with unlimited resources (infinite 
carrying capacity), the population grows exponentially. 
The Malthus model is an example that generates this be- 
havior. However, for finite carrying capacity, the growth 
of individual organisms Q , tumours [§[ and other biologi- 
cal systems @ are well fitted by sigmoid curves [J-Q that 
can be obtained from the Gompertz or Verhulst model, 
for instance. 
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The generalization of the well known growth models, 
by the von Foerster et al. [Iol | or Richards [llj ] models, 
introduces additional parameters to the former models, 
making them more complete and accurate to fit exper- 
imental data. In this attempt to find a suitable com- 
plete growth model, the one-parameter generalizations 
of the lo gar ithm and exponential functions play a central 
role [12h16| |. They allow us to easily retrieve particu- 
lar cases without needing to calculate limits and permit 
convenient algebraic tricks to handle the expressions. 

From the analytical solutions obtained from the one- 
species growth models that we address, we call attention 
to the following items. For each model, the population 
evolution is proportional to the inverse of the growth 
rate. This gives rise to a dimensionless characteristic 
time and can be defined as the system independent vari- 
able r. Furthermore, we define a scaling function that 
depends on the solution of each model. The dependent 
variable y is the ratio between the scaling function, at a 
given time r, with respect to its initial value. We show 
that using these variables, the models are independent 
from the initial conditions and parameters demonstrat- 
ing a universal Malthusian behavior. This universality 
occurs for a finite carrying capacity models and even in 
the Tsoularis- Wallace model, where no closed analytical 
solution is found. 

Finally, we stress that when a constant effort rate is 
considered, the universality is preserved and the steady 
state (asymptotic) solution can be interpreted as an order 
parameter. A transition between extinction and survival 
phases are separated by a critical value, which depends 
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on the effort rate. This transition does not occur for 
the Gompertz model. This picture is not altered when 
time dependent effort rates is considered in the Richards' 
model, nor when a time-dependent growth rate is consid- 
ered. Although a closed analytical solution is not found, 
we are able to show the universality of the Tsoularis- 
Wallace model with constant effort rate. 

This manuscript is structured as follows. In Sec. [ill 
we briefly review the one-parameter generalization of the 
logarithm and exponential functions and present some of 
its properties used along the study. In Sec. Mil we intro- 
duce the main one-species population dynamics (growth) 
models. We define the scaling function and show that all 
these models can be written as the universal Malthus (ex- 
ponential) model. In Sec. IIV1 we consider the insertion 
or removal of individuals through a constant effort rate. 
We show that this quantity does not affect the univer- 
sality of the models. Moreover, it induces an extinction- 
survival transition at a well determined value. Next, we 
show that time dependent effort rate and time depen- 
dent growth rate do not affect universality. In Sec. |Vl we 
present a formula to obtain the scaling functions and our 
final remarks. 



III. GROWTH MODELS 

The one-species growth models can be characterized 
by the saturation function: 

G(N) = dh ™ , (3) 
dr 

which is the per capita growth rate. Time is measured 
as the inverse of the growth rate i.e., r = nt. 

Below, we briefly introduce the basic one-species popu- 
lation dynamics (growth) models. These models fall into 
two categories: one with infinite (described in terms of 
the number of individuals N) and the other with infinite 
(p = NJ K) carrying capacity. 

A. Infinite Carrying Capacity 

We start presenting the Malthus model and a non- 
linear gen eralization given by the von Foerster's et al. 
model MIS- The y 

are described in terms of the num- 
ber of individuals N since they have infinite carrying ca- 
pacities. 



II. GENERALIZED LOGARITHMIC AND 
EXPONENTIAL FUNCTIONS 



In the following, we introduce the one-parameter gen- 
eralization of the logarithmic exponential function and 
present some of their main used properties. 

The q-logarithm function is defined as: 



ln^(x) 



1 



dt 



(1) 



This one-parameter generalization of the natural loga- 
rithm function, which is retrieved for q — > 0, has been 
introduced in the context of non-extensive statistical me- 
chanics [H, [f| and is defined as the value of the area un- 
derneath the non-symmetric hyperbole, fq(t) = l/t 1 ^^, 
in the interval t G [l,x] [Hj]. Note that in Eq. (fTJ), 
ln.q(x) is not "logarithm x in the base <f. For q < 0, 
liiq(oo) = — 1/q; for q > 0, lrXq(O) = — 1/q; for all q, 
lnq(I) = 0; lnq(x _1 ) = — ln_g(x); dlriq(x) / dx — x^ 1 . 

The inverse of the q-logarithm function is the q- 
exponential function 



Sg(x) 



liny_, 5 (l + q x) ■ 




if qx > — I 
otherwise 



(2) 



so that e q (0) = 1, for all q and [e q 



eq/ a (ax), where 



a is a constant. 



1. Malthus Model 



For the Malthus model G{N) = 1 
dN{r) 



dr 



N(t) 



(4) 



the result is a populational exponential growth: N(r) 
Noe T . The scaling function is: 



writing 



s (t)=N(t) 



, \ s (r) 
y(T) = ~ 7~r = e 



»o(0) 



(5) 



(6) 



one obtains the universal equation, that is independent 
from the parameter k (growth rate) and initial condition 
N . 



2. von Foerster et al. Model 

In the von Foerster's et al. model OH : G a (N) = N a 
and: 



dN(r) 
dr 



N(t) 



l + Q 



(7) 



where a is the generalization parameter. Its solution is: 

1 



N(r) 



[a(T- T )]V a ' 



(8) 
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where T is a dimensionless time at which the population 
size diverges. As a —> 0, one retrieves the Malthus model: 
N(t) = Noe T , with No = e T . For human population, 
a fs 2 and a doomsday is predicted to occur on a Friday, 
November 13^, 2026 [l(|. In terms of the one-parameter 
of the generalized exponential function, Eq. ([5]) can be 
written as [l7| : 



N(t) 



(9) 



No e Q (-iV «r) 

where N = N(0) = l/(aT) 1 /". Figure ©(a) depicts 
N(t) for three different set of parameters. 

Consider Eq. M as ln_ Q [iV(T)/7V ] = N°t, using 
Eq. Q] one obtains: N^ a — N~ a (r) = —ar, after re- 
grouping it, one writes ln_ Q [iV(T)] — ln_ Q [iVo] = r or 



,ln_ Q [Ar(r)] / e ln_ Q [JV ] 



e T . In this way, we are able to 



define the scaling function as: 



s a (r) 



Jn- a [N(r)] 



(10) 



so that one retrieves the Malthus scaling function (Eq. [3]) 
for a — 0. Writing y = s Q (r)/s Q (0) = e T , one sees 
that it becomes independent from the parameters and 
from the initial condition, presenting universality. This 
universality is depicted in Figure ©(b). 



B. Finite Carrying Capacity 

The environment with limited resources is taken into 
account in the dependent variable p(r) = N{t)/K. Opti- 
mum environment exploration is achieved when p{r) = 1. 
Given the initial condition t = tq — and po = p(0). 
The steady state solutions, according to their stability, 



[dp(T)/di 



= p*G{p*)] are either species extinction 



p " = or survival p*, obtained as the roots of G(p*) = 0. 

In the following, we present the Gompertz and Verhulst 
models and their generalization, the Richards' model [ill ]. 
One way to include the finite carrying capacity K to the 
von Foerster et al. model is to replace the number of 
individuals N by ]n(N/K) in the saturation function, 
this leads to the hyper-Gompertz model [III [l!| . These 
models have closed analytical solutions and we can col- 
lapse the solutions on a single exponential equation. The 
Tsoularis- Wallace model [20| generalizes all these models, 
but it does not have a general analytical closed solution. 
However, we are able to write a scaling function and show 
that its solution is universal. 



1. Gompertz Model 

The Gompertz model, Go(p) = — hip has the solution: 

p(r) — e( lnp °) e . Writing lnp(r)/lnpo = e ~ T , one ob- 
tains the scaling function: 

s (r) = {Hpir)}}- 1 . (11) 

Defining y = sq(t) / sq(0) = e T , one retrieves the univer- 
sal Malthusian curve. 
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Figure 1. (a): Evolution of the population according to the 
von Foesters model (Eq. [9| for three different set of parame- 
ters: (I) iVo = 10; a = 1.0; k = 0.05, note that the divergence 
occurs for t = 2.0 (II) 7V = 50; a = 0.0 (Malthus); k = 0.4, 
(III) iVo — 30; ex — — 1.0; k = -10.0. (b): The universal 
curve y = s a (r) / 'sf(O) = e T , with s a given by Eq. [TUJ Inset: 
Zoom in of the specified region. 



2. Verhulst Model 

The Verhulst model G\(p) = 1 — p has the following 
solution: 



p(r) 



1 



1-{1-P^)t 



(12) 



Since r = nt, the model is independent from the pa- 
rameter k but still depends on the initial condition pq. 
Frequently, authors assign a partial universal behavior to 
the Verhulst equation: it is independent only on the pa- 
rameter k, but not on initial condition pq. Nevertheless, 
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calling the scaling function as 

Si (t) = [1- P -\t)}- 1 



(13) 



one can write the ratio of the scaling function with re- 
spect to its initial values 



y 



*i(t) = [l-p-Hr)] 
Sl (0) 



Po 1 }- 1 



(14) 



With this procedure, one gets rid of the dependence on 
the initial condition and retrieves full universality. E. W. 
Montroll used this transformation in Ref [2lJ, p. 4634] but 
he has not called the attention to its universal aspect. 



3. Richards' Model 

Now consider the Richard's model G q (p) = — (p' — 
[Hi EE HH HI]) which can conveniently be written 
in terms q- logarithmic function (Eq. [IJ: 



dlnp(r) 



In," [p(r)] 



(15) 



In Ref. [23[, the authors assume that the replication 
rate of a cell is regulated by a competition between the 
cell impetus to proliferate and an inhibition from the 
other cells through diffusive growth factors. These as- 
sumptions result in a differential equation for the growth 
of the cellular system identical to EairBI where the pa- 
rameter q is related to the range of interaction between 
cells and the fractal dimension where these cells grow [l5[ . 

Notice the analogy of Eq. [13] to the Gompertz model, 
where we have changed the natural logarithmic function 
in the saturation function by the generalized one. When 
discretized, this equation leads to the generalized logistic 
map [l6j |. For q = 0, one retrieves the Gompertz model 
and for q = 1, the Verhulst one. The solution of Eq. (fTS"]) 
can be written in terms of the generalized exponential 
and logarithm functions as 



p(r) 



1 



e 9 [ ln 9 (PO X ) e ~ 



(16) 



The asymptotic limit (r — > oo) of Eq. fTf)]) is p* — 
p(oo) = 1, regardless the choice of q. 

Let us now analyze the universality of the Richards' 
model. From Eq. (|TB| and using the properties 
of the generalized functions, one can write p(r) = 
e_g[— In, (pq 1 ) e~ r ] = e_q[ln_q (p ) e~ T ]. Nevertheless, 
taking the —q- logarithm of the preceding equation, one 
obtains ln g [p(r)] / ln 9 [p ] = e _T . Writing the scaling 
function as 



(17) 



one generalizes the Gompertz (q = 0, Eq. [TT|) and 
Verhulst (q = 1, Eq. [T3"|) scaling functions. Defining 
y = Sq(r)/ Sq(0) = e T , the Richards model presents a 
universal behavior, since the model is independent from 
the parameters (q, k and K) and from the initial condi- 
tion (p ). 



4- The Hyper- Gompertz Model 

The hyper-Gompertz model is obtained from the von 
Foerster et al. model, replacing N a by [—ln(N/K)]~ t 
and calling p = NIK, the saturation function becomes 
G 7 (p) = [-lnp]T >6], so that [ll[igj]: 



dlnp(r) 
aW 



{-ln[p(r)]r, 



(18) 



which solution is: 

lnp(r) = -[( 7 - 1)t - (-1)*" (lnpo) 1 - 7 ] 1 ^-^ , (19) 
which can be written as: 

[- lnp(r)] 1 -^ - [- lnpo] 1 " 7 = (7 - 1)t , (20) 
calling the scaling function as: 

fl7 (r) = e[- to ^ 1 " 7 /<T- 1 ), (21) 
the ratio y = s 7 (r)/s 7 (0) = e T is universal. 

5. Tsoularis-Wallace Model 

So far, all the presented models have closed analytical 
solutions. Let us now unify these models, which presents 
an analytical solution, but not a general closed form. In 
terms of the g-logarithm function, the Tsoularis-Wallace 
model, G aAn (p) =p a [-ln 9 -(p)F, is M- 



dlnp 
dr 



.p a (T){-ln q [p(r)}y . 



(22) 



The Richards' model [Eq. [TS] is retrieved for a = and 
7 = 1. For a = q = 0, one retrieves the hyper-Gompertz 
model [Eq. [15]. Writing p = N/K, 7 = and rescaling 
the growth rate, one retrieves von Foerster et al. model 
[Eq.0. 

The solution p(r) of Eq. [22] is the root of: 
Bp« (r ) f-| , 1 - l)-B pt (-- 1 - 7 ) = <f-"V , (23) 



where B x (a,b) — J Q t a 1 (1 — t) b 1 dt is the incomplete 
beta function. Writing the scaling function as: 



(24) 



for y — s a t q , 7 (r)/s aj § , 7 (0) = e r , one obtain the univer- 
sality of the Tsoularis-Wallace model. This model il- 
lustrates that to find the scaling function, one does not 
imperatively has to know the analytical closed solution 
for p(r) . 
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IV. EXTINCTION AND SURVIVAL PHASES: 
EFFORT RATE 

The effort rate e quantifies insertion and removai of in- 
dividuais in a population. Here, we investigate the effect 
of e in the universality of Richards (closed analytical solu- 
tion) and Tousalis- Wallace (analytical solution) models. 
Further, we point out the transition between the extinc- 
tion and survival phases at a determined critical value. 
Although we do not consider the stochastic models, the 
time dependent effort rate can be considered as a ran- 
dom variable e(i). In this case, one is able to deal with 
additive noise. To treat multiplicative noise, one must 
consider a time dependent growth rate n{t). 




A. Constant Effort Rate 

With constant effort rate, one can deal with the 
transition between the extinction and survival phases. 
The steady-state solution represents the order parameter 
since it vanishes in the extinction phase and, at a critical 
point, it gives rise to a survival phase. This transition is 
present in the Richards and Tsoularis- Wallace models. 



1. Richards- Schaefer's Model 

We refer to the Richards' model (TTSl) with a constant 
effort rate as the Richards-Schaefer's model: 



d\np(r) 



•lng[p(r)] +e , 



where e = e/n. The solution of the Eq. (|2"5]) is: 



P (t) 



°S { ln <? 



y 



(25) 



(26) 



where po = p(0) is the initial condition. In Fig. @, the 
different behaviors of p(r) for different set of parameters 
are depicted. For e = in Eq. (|2l)l) . one retrieves the 
Richards' model Eq. [I])] 

The asymptotic behavior of Eq. (12T))) is: 



p* =p{oo) = e,(e) 



(27) 



so that for e = 0, one retrieves p* = 1, as expected 
from Eq. [26l It is interesting to point out that Eq. ((27)) 
vanishes for qe < —1, representing species extinction (see 
Fig. |3J). Species survival occurs for e > e^ c \ where the 
critical value is: 



q 



(28) 



One thinks of p* as an order parameter that describes two 
ecological stable phases: species extinction and survival. 
These phases are separated by the critical value e c . This 
transition is supressed for q — 0, the Gompertz model. 
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Figure 2. (a): Evolution of the population according to the 
Richards-Schaefer's model (Eq. l26[) as a function of r for three 
different set of parameters: (I) po = 0.1; q = 2.0; e = —0.1, 
(II) po = 0.6; q = 0.5; e = -0.75,(111) p = 0.8; q = 1.0; e = 
— 2.0.(b): the universal curve y = Sq(r)/sq(0) — e T , with Sq 
given by Eq. 1291 Inset: Zoom in of the specified region. 



To show the universality of this model, write 
y(r) — p(T)/eq(e), so that Eq. ([2"rJ| becomes 
y(r) = e_g{{ln_g(y )exp{[-(l + qe)r}}, where y = 
?/(0) and we have used the properties: l/eq{x) — 
e-q(-x) and — ln^x -1 ) = ln_g(x). Apply- 
ing — q-logarithm on y(r) one has: ln_g[?/(T)] = 
ln_g(y )eg(e) exp{[— (1 + qe)r}. Thus Eq. [56] can be 



written as ln_ ? [p(r)/e ? (e)]/ ln_^[p /e,(e)] 
The scaling function is given by: 



ln-9 [P( r )/ e q 



e q (e) 



(29) 



So that e = => e^(0) = 1 and one retrieves Eq. [TTJ 
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B. Time Dependent Effort Rate 
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Figure 3. Asymptotic behavior of the Richards-Schaefer's 
Model solution, given by Eq. ([27]), as a function of e for sev- 
eral values of q. The species extinction (p* = 0) happen only 
when q > 0. The critical e values, i.e. the e value which 
extinction take place, is given by Eq|2"51 



Calling y = Sq(r)/ Sg(0) — e T , one retrieves the universal 
Malthus model. The introduction of e^(e) in the denom- 
inator of the scaling function, which does not affect the 
result, is due the necessity of it to be compatible with 
the time dependent effort rate. 



2. Tsoularis-Wallace-Schaefer Model 

Let us consider an effort rate in the Tsoularis- Wallace 
model, and call it the Tsoularis-Wallace-Schaefer model: 



dh\p(r) 
~d~r 



= p<*(T){-ln s \p(T)]y+. 



which solution p(r) is the root of 
«* T > dx 



x{x a [- hxq{x)]i + e} 
The scaling function is defined as 

/• p(T) dx 



"a, q,7 



, e (r) = exp 



x{x a [- lHg(x)]T + e} 



(30) 



(31) 



(32) 



so that a universal behavior is found for 



;(0) = e T , even so we have not been able 



to calculate the integral to find explicitly p(r). 

The steady-state solution is obtained considering 
dp/dr = in Eq. HB so that p*{p* a [- ln,-(p*)] 7 + e} = 
0, leading to p* = 0, which represents the extinction 
phase and p* a [~hiq(p*)] 1 — — e, which represents the 
survival phase. The critical value separationg the ex- 
tinction and survival phase, is given by the root of: 
(p*)«/7 _ ( p *)«/7+§ = q(-e) 1 ^. 



It is interesting to consider time dependent effort rate 
once it can be considered as a random variable and noise 
can be treated. On one hand, if the growth rate is con- 
stant, the noise is additive to the model. On the other 
hand, if the growth rate is time dependent and a random 
variable, one can consider the effort rate to vanish and 
one has multiplicative noise in the system. 

Although we do not address the stochastic models, we 
present the analytical solutions of these models and show 
their universalities. 



1. Constant Growth Rate 

Consider now a time dependent effort rate e(i) in the 
Richards-Schaefer's model [Eq. (l2"5j)]: 



dlnp(r) 
~d~T 

The solution of Eq. 
p(r) = 



= -faqP(r) +e(r) . 



is 



where 



e( T ) = I f dT'eir' 

T JO 



(33) 



(34) 



(35) 



is the mean value of e(r) up to time r. For a constant 
effort rate e(r) = e in Eq. [MJ one retrieves the Richards- 
Schaefer's model and its solution [Eq. (J2BJ], 

The steady state solution (r — > oo) of Eq. IMlis: 



p* =p(oo) = eg(e) 



(36) 



where e = e(oo) is the true mean value of e(r). Species 
extinction occurs for qe < —1. 

Write y(r) = p(r)/e $ [e(r)], so that Eq. 



becomes y(r) = e_ 5 {{ln_q(y )/e§[e(0)]} 



<r)] 



exp{ — [1 + ge(r)]r}}, where yo — y(0) and we have 
used: l/eq(x) = e_^(— x) and — ln^x -1 ) = ln_g(x). 
Applying — g-logarithm on y(r) one has: ln_g[y(r)] = 
ln_g(y o )eg[e(r)]/e 5 [e(0)] exp {-[1 + ge(r)]r}, justifing 
the definition of the scaling function as 



ln_g 



e f [e(r)] 



eMr)] 



(37) 



Taking e(r) = e(0) = e, one retrieves Eq. [23 Notice 
that it becomes clear the reason we have introduced the 
factor l/eq(e) in the definition of the potential growth 
on the Richard-Schacfcr model with constant effort rate 
Eq.EH 
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Figure 4. The considered growth models from sections II II I and IIV Al with their respective quantities: Saturation function G 
and scaling function s. The Malthusian universal behavior (independence of parameters and initial conditions) is obtained in 
all models by the ratio of the scaling function in relation to its initial value: y = s(t)/s(0) = e T . 



Defining y — Sq, e ( r )( T )/sg ie ( T )(0) = e r , one obtains the 
universal Malthusian equation, independent of parame- 
ters and initial conditions. 

If e(t) is a random variable, then one has the additive 
stochastic growth equation. In this case, if its mean value 
vanishes e(r) = and e(Ti)e(T 2 ) = ct 2 <5(t 2 — t\) (Gaus- 
sian process), then the probability density function of 
v = lap satisfies the Fokker-Planck equation: d T P(v) = 
d v [P(v)\n <j (v)} + (a 2 /2)d^[P(v)} Correlated 
and Levy like noise have also been addessed [26l - l28j . 



2. Time Dependent Growth Rate 

Consider now time dependence in both growth and ef- 
fort rates. One has: 



dHp(t)} 
dt 



-re(t) ln§[p(t)]+e(t). 



(38) 



Notice that here we use t instead of r as the independent 
variable. The soultion of Eq. [35] is given by: 



dt'l{t')K{t')+p-' 



-1/9 



where: 

I(t) = e U? dt' K (t')+g ft dt'e(t')] 
so that 1(0) = 1. 



(39) 



(40) 



i(t) P -"(t) - i(o) Po q = / dt'i(t'Mt') (4i) 



so that 



e i{t)p- q (t) 



= /* dt'i{t') K {t') 



(42) 



Calling f = L dt' I(t')ri(t') and the scaling function as: 

M0,5,e-(t)(t) = e f « p " 4 W. (43) 

The ratio y = s re (t) )$) g( t )(t)/s re ( t ) )$ig ( t )(0) = e T is univer- 
sal. 

If e(t) = 0, one can now consider the time dependend 
growing rate as n(t) = cir,(t) + ai7i(i), where ao(t) a de- 
terminist growth and 71 (t) ma y b e considered as a mul 
tiplicative stochastic noise 



ma y b e 
[29.'30]. 



V. CONCLUSION 

We show that measuring time r = tj 7c, where k is the 
growth rate, we are able to write a general expression to 
obatain the scaling function: 



S (r)=e/»*/ G ( e "', 



(44) 



where, v = In AT and G(N) is the saturation function. 
Using the scaling function, all the considered models are 



8 



written as the universal Malthus (exponential) model. 
We have shown the universal properties of the Tsoularis- 
Wallace-Schaerfer model, which is a very general con- 
stant coefficient model, with no closed analytical solu- 
tion. Although we have not derived all the possible par- 
ticular cases from this model, we can deduce that they are 
all universal. If one includes the effort rate, a transition 
from species extinction to the survival is well determined. 
This effort rate may represent the mean field approxima- 
tion of the interaction of other species. For this reason, 
we believe in the universality of multi-species models. 
For time dependent coefficients, the most general model 
we have addressed and solved is the Richards-Schaeffer 
model, which is also universal. Since one can consider 
stochasticity in the time dependent coefficient models, ei- 
ther with additive and multiplicative noise, we conjecture 



that the stochastic models are also universal. The uni- 
versality is useful when working with experimental data 
as the model that best fits the data must correspond to 
a straight line in a data collapsed semi-log graph. 
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